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Abstract 

We introduce a new method to treat Majorana fermions on the GRACE system 
which has been developed for the automatic computation of the matrix elements for 
the processes of the standard model. In the GRACE system, we already have such 
particles as Dirac fermions, gauge bosons and scalar bosons within the standard 
model. On the other hand, in the SUSY models there are Majorana fermions. In 
the first instance, we have constructed a system for the automatic computation of 
cross-sections for the processes of the SUSY QED. The system has also been applied 
to another model including Majorana fermions, the minimal SUSY standard model 
(MSSM), by the re-definition of the model file. 



1 Introduction 

It has been a promising hypothesis that there exists a symmetry called supersymmetry 
(SUSY) between bosons and fermions at the unification-energy scale. It, however, is 
a broken symmetry at the electroweak-energy scale. The relic of SUSY is expected to 
remain as a rich spectrum of SUSY particles, partners of usual matter fermions, gauge 
bosons and Higgs scalars, named sfermions, gauginos and higgsinos, respectively M. 

The quest of these SUSY particles has already been one of the most important pursuits 
to the present high-energy physics || . Although such particles have not yet been discov- 
ered, masses of them are expected to be O(10 2 ) GeV ||. In order to obtain signatures 
of the SUSY-particle production, electron-positron colliding experiments are preferable 
because the electroweak interactions are clean and well-known. Thus we hope SUSY 
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particles will be found out at future TeV-region (sub-TeV region) e e + -colliders such as 
CLIC, NLC and JLC [§. 

For the simulations of the experiments, we have to calculate the cross-sections for the 
processes with the final 3-body or more. We have already known within the standard 
model that the calculation of the helicity amplitudes is more advantageous to such a 
case than that of the traces for the gamma matrices with REDUCE P, H|. The program 
package CHANEL is one of the utilities for the numerical calculation of the helicity 
amplitudes. 

It, however, is also hard work to construct a program with many subroutine calls of 
CHANEL by hand. Thus we need a more convenient way to carry out such a work. The 
GRACE system ||, which automatically generates the source code for CHANEL, is one of 
the solutions. The system also includes the interface and the library of CHANEL, and the 
multi-dimensional integration and event-generation package BASES/SPRING v5.1 ||. 

In the SUSY models, there exist Majorana fermions as the neutral gauginos and hig- 
gsinos, which become the mixed states called neutralinos. Since anti-particles of Majorana 
fermions are themselves, there exists so-called 'Major ana-flip', the transition between par- 
ticle and anti-particle. This has been the most important problem which we should solve 
when we realize the automatic system for computation of the SUSY processes. 



In a recent work [10|, [Tl|] , we developed an algorithm to treat Majorana fermions in 
CHANEL. In the standard model, we already have such particles as Dirac fermions, gauge 
bosons and scalar bosons in the GRACE system. Thus we can construct an automatic 
system for the computation of the SUSY processes by the algorithm above in the GRACE 
system. In this work, we present the check list of the system at this time, and one of the 
results. 



2 Majorana fermions into new GRACE 

In Fig. 1, we present the system flow of GRACE [|12|]. The GRACE system has become 



more flexible for the extension in the new version called 'grc' |[L3| , which includes a new 
graph-generation package. With this package, every graphs can be generated based on 
a user-defined model. It is necessary for us to make the interface and the library of 
CHANEL and the model file for including the SUSY particles. 

The method of computation in the program package CHANEL is as follows: 

1. To divide a helicity amplitude into vertex amplitudes. 

2. To calculate each vertex amplitude numerically as a complex number. 

3. To reconstruct of them with the polarization sum, and calculate the helicity ampli- 
tudes numerically. 

The merit of this method is that the extension of the package is easy, and that each vertex 
can be defined only by the type of concerned particles. 

Here we propose an algorithm [|K], for the implementation of the embedding Ma- 



jorana fermions in CHANEL as follows: 
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Fig. 1. GRACE system flow 
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• policy 

1. To calculate a helicity amplitude numerically. 

2. To replace each propagator by wave functions or polarization vectors, and 
calculate vertex amplitudes. 

3. Not to move charge-conjugation matrices. 

• method 

1. To choose a direction on a fermion line. 

2. To put wave functions, vertices and propagators along the direction in such a 
way: 

i) To take the transpose for the reverse direction of fermions 

ii) To use the propagator with the charge- conjugation matrix for 
the Majorana-flipped one. 

As a result, the kinds of the Dirac-Majorana-scalar vertices are limited to four types: 

(1) UTU 

(2) U T TU T 

(3) Z7c T r T F T 

(4) XFV^C^U 

where U's denote wave functions symbolically without their indices, and C is the charge- 
conjugation matrix. The symbol T stands for the scalar vertex such as 

r = *.iz2 + JlB .i±2 . 

The vertices (2)~(3) are related to the vertex (1) which is the same as the Dirac-Dirac- 
scalar vertex in the subroutine of CHANEL. Thus we can build three new subroutines for 
the added vertices. We have performed the installation of the subroutines above with the 
interface on the new GRACE system. 



3 Numerical results 

At the start for the check of our system, we have written the model file of the SUSY 
QED. In this case, there is only one Majorana fermion, photino. Next we have extended 
the model file and the definition file of couplings for the MSSM. The tests have been 
performed by the exact calculations with the two methods, our system and REDUCE. 
In Table I, the tested processes are shown as a list. The references in the table (without 



IIJ) are not the results of the tests, but for help. 

Here we present the results for the single-selectron production within the SUSY QED 
process at the energy y/s = 190 GeV. The masses of the concerned particles are = 50 
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Process 



Number of diagrams Comment 



SUSY QED 



Check Reference 



e e 



e R e R 
e R e L 



2 

2 
2 



Major ana-flip 
in internal lines 



OK 
OK 
OK 



El 



e e n 



e R e R 



Including pair 
annihilation 



OK 
OK 



El 



El 



e e" 1 



e R e L 
e R e L 



Values are 
equal 



OK 
OK 



El 
3 



e e n 



77 



F-B symmetric 



OK 



m 



e e" 1 



777 



12 



Final 3-body 



OK 



e e n 



%7e 



12 



Including every 
elements for tests 



OK 



MSSM 



e e 



e T e 



4 Majorana fermions OK 



e e 1 



Xi Xi 



OK 



Table I. The list of the tested processes. 



GeV, Mg H = 100 GeV and Mg L = 130 GeV in the calculation. This is the case that the 
pair-production processes occur for both selectrons at the JLC-I energy, but they do not 
at the LEP-II energy. The Feynman diagrams for this process, which are drawn by the 
program package 'gracefig' [Tj| in the new GRACE, are shown in Fig. 2. 



In Fig. 3, we show the angular distribution of the outgoing positron in the process 
e~e + — > e^7e + . Here we use BASES for the calculation from the REDUCE output. The 
result is in beautiful agreement with the value that is obtained by GRACE at each bin 
of the histogram. Since the two diagrams with the one-photon exchange dominate in this 
case, there is a steep peak in the direction of the initial positron. In such a case, the 
equivalent-photon approximation (EPA) works well JT7|| . 

In Fig. 4, we show the T$ distribution of the selectron in the process e~e + — > e^ ; ye + . 
The quantity T s is defined in Ref. |T7] 



as 



P T 2 + M gR 2 

Ts - m 2 ' 



where Pt denotes the transverse momentum of the selectron. We show also the Pt 
distribution of the selectron in the process e~e + — > e^'je + in Fig. 5. Here we calculate 
the two dominant diagrams for comparison with the result from EPA. 



4 Summary 

We introduce a new method to treat Majorana fermions on the GRACE system for the 
automatic computation of the matrix elements for the processes of the SUSY models. In 
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the first instance, we have constructed the system for the processes of the SUSY QED 
because we should test our algorithm for the simplest case. The numerical results convince 
us that our algorithm is correct. 

It is remarkable that our system is also applicable to another model including Ma- 
jorana fermions (e.g. the MSSM) once the definition of the model file is given. We 
have calculated the processes e~e + — > 777 and e~e + — > e^f)e + within the SUSY QED. 
We should calculate the single-photon event from e~e + — > X1X1I EH; an d the resultant 
single-positron (electron) event from the single-selectron production e~e + — > e T Xi e± [0 
as soon as possible. It should be emphasized that the GRACE system including SUSY 
particles is the powerful tool for this purpose. 
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